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SCATTERING  PARAMETER  TRANSIENT 
ANALYSIS  OF  TRANSMISSION  LINES 
LOADED  WITH  NONLINEAR  TERMINATIONS 


Abstract 

This  work  presents  a  new  approach  for  the  time-domain  simulation  of  tran¬ 
sients  on  a  dispersive  and  lossy  transmission  line  terminated  with  active  devices. 
The  method  combines  the  scattering  matrix  of  an  arbitrary  line  and  the  nonlinear 
causal  impedance  fuiictions  at  the  load  to  derive  a  closed-form  algorithm  for  the  sig¬ 
nal  at  the  near  and  far  ends. 

The  problems  of  line  losses,  dispersion  and  nonlinearities  are  first  investigated. 
A  tune-domain  formulation  is  then  proposed  using  the  scattering  matrix  representa¬ 
tion.  The  algorithm  assumes  that  dispersion  and  loss  models  for  the  transmission 
lines  are  available  and  that  the  frequency  dependence  is  known.  Large-signal 
equivalent  circuits  for  the  terminations  are  assumed  to  be  given.  Experimental  and 
computer  simulated  results  are  compared  for  the  lossless  dispersionless  case,  and  the 
effects  of  losses  and  dispersion  are  predicted. 
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Figure  1.  Transmission  line  with  nonlinear  terminations 

and  source  generator . 

Figure  2.  Transmission  line  and  frequency  domain  flow-graph 

representation  using  scattering  parameters . 

Figure  3.  Time-domain  circuit  and  flow-graph  representation 
of  arbitrary  transmission  line  (TL)  with  nonlinear 
terminations  at  time  t.  The  sign  *  indicates  a 
convolution  between  time-domain  scattering  parameters 
and  the  independent  voltage  waves  ai  and  a2  as  per 
Equations  (14)  and  (15)  . 


Figure  4.  Equivalent  flow-graph  representation  of  arbitrary 
transmission  line  with  nonlinear  time-varying 
terminations  at  time  t.  The  history  of  the  network 
is  included  via  the  memory  variables  M^(t)  and  M2(t). 
Note  that  this  representation  involves  only 
multiplications  in  the  time  domain . 


Figure  5.  Linear  capacitor  and  equivalent  time-domain 

Backward-Euler  representation  . 

Figure  6.  Geometrical  and  circuit  interpretation  of  the 

Newton-Raphson  algorithm  at  the  ktl*  iteration  .  .  .  . 


Figure  7.  Comparison  of  theoretical  (plots)  and  experimental 
(photographs)  simulations  for  a  stripline  structure 
terminated  with  ALS  inverters  at  near  (left)  and 
far  (right)  ends.  The  line  length  is  50  inches. 

The  low-frequency  characteristics  are  ZQ  -  730  and 
v  »  0.169/ns.  Pulse  width  is  160  ns;  rise  and  fall 
times  are  4  ns.....  . 


Figure  8.  Comparisons  of  simulated  waveforms  for  a  lossless 
(solid  line)  and  lossy  (dashed  line)  microstrip 
line  with  500  at  the  near  end  and  open  at  the  far 
end.  Line  characteristics:  length  -  25  inches, 
low-frequency  characteristic  impedance  :  Zn  »  600 
with  loss  characteristics  :  R  “  500/(ra-  /2irGHz) , 

G  *  0.02mhos/m.  Pulse  characteristics  :  width  *  20ns, 
rise  and  fall  times  :  4ns.  In  both  cases,  the 
dispersion  model  is  that  of  [3]  . 
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I.  INTRODUCTION 

In  today’s  many  applications  of  integrated  circuits  and  printed  circuit  boards, 
transmission  lines  and  interconnections  play  an  instrumental  role  at  virtually 
every  level  of  integration.  With  the  design  of  fast  devices  having  switching  times 
in  the  picosecond  range,  transmitting  data  at  high  megabaud  rates  has  become  very 
commonplace  in  modern  digital  computers  and  switching  networks  used  for 
telecommunication.  Signal  delays  and  rise  times  are  more  and  more  limited  by 
interconnection  lengths  rather  than  by  device  speed  and  represent  a  potential  obs¬ 
tacle  to  the  ultimate  scaling  on  VLSI  technolgy.  In  recent  years,  modeling  inter¬ 
connections  has  become  a  major  focus  of  interest  in  the  implementation  of  digital 
and  microwave  circuits.  Shorter  rise  and  fall  times  as  well  as  higher  frequency 
signals  have  compelled  most  transmission  lines  to  operate  within  ranges  where 
dispersion  is  no  longer  negligible.  Skin  effect  and  losses  contribute  to  signal  corrup¬ 
tion  leading  to  waveform  attenuation  and  a  reduction  of  the  effective  dimensions. 
In  wafer  scale  integration,  these  losses  can  become  very  significant  and  may  lead  to 
an  RC  type  behavior  of  the  lines.  Finally,  in  the  case  of  multiconductor  lines, 
cross-coupling  between  neighboring  lines  may  increase  the  level  of  distortion  in 
excited  lines  which  can  initiate  false  signals  in  non-excited  lines. 

The  implementation  of  a  high-density  compatible  packaging  scheme  is  essen¬ 
tial  for  the  design  of  high-speed  digital  systems  such  as  gallium  arsenide  integrated 
circuits,  for  microwave  or  digital  applications,  printed  circuit  boards,  chip  carriers, 
and  modeling  of  these  networks  represents  the  first  step  toward  implementing  reli¬ 
able  design  guidelines.  A  complete  CAD  tix.il  for  studying  these  effects  would 
require  a  frequency  domain  characterization  of  the  transmission  line  with  higher- 
order  modes  included  to  account  for  dispersion.  Numerous  authors  have  investi¬ 
gated  ihe  properties  of  microstnp  lines  at  high  frequencies  and  derived  expressions 
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relating  the  propagation  characteristics  to  frequency  IlHlO].  Full-wave  and 
simplified  models  have  been  proposed  to  describe  these  effects  and  derive  the  fre¬ 
quency  dependence  of  the  characteristic  impedance  and  the  propagation  constant. 
Other  geometries  such  as  stripline,  buried  microstrip  or  coplanar  have  thus  far 
received  less  attention  but  obey  the  same  restrictions  imposed  on  the  electrical  per¬ 
formance  of  microstrip  at  microwave  frequencies. 

Of  equal  importance  is  the  analysis  of  a  high-speed  or  high-frequency  signal 
propagating  on  a  dispersive  and  lossy  transmission  line.  Such  an  analysis  requires  a 
complete  and  accurate  frequency  characterization  of  the  structure  of  interest  and, 
for  practicality,  must  implement  the  nonlinear  and  time-changing  behaviors  of  the 
terminations,  which  are  transistors,  logic  gates  or  other  types  of  active  devices. 
Several  investigators  have  attempted  to  set  up  analytical  models  describing  wave 
propagation  in  such  systems.  Solutions  for  lossless  lines  with  arbitrary  termina¬ 
tions  were  obtained  by  Mohammadian  et  ed.  [12]  using  a  forward  and  backward 
wave  approach.  Veghte  and  Balanis  [  1 3]  have  analyzed  the  distortion  of  a  pulse  due 
to  dispersion  along  a  microstrip  transmission  line.  Djordjevic  et  ul.  [  1 1  ]  used  a 
complementary  network  approach  to  simulate  the  time-domain  transient  on  a  mul- 
t [conductor  array  with  nonlinear  terminations.  However,  exact  closed-form  solu¬ 
tions  that  combine  losses  and  nonlinear  terminations  have  not  been  made  available 
as  yet. 

In  this  study,  a  combined  frequency-domain  time-domam  approach  is  used  to 
formulate  the  propagation  equations  on  a  dispersive  and  lossy  line  with  nonlinear 
behavior  at  the  terminations.  The  method  uses  the  scattering  parameters  of  the 
transmission  line  to  derive  a  closed-form  algorithm  for  the  propagation  ot  an  arbi¬ 
trary  signal  in  the  time  domain.  In  the  lossless  case,  the  solution  reduces  to  very 
simple  expressions  which  greatly  increase  computational  efficiency. 


II.  FORMULATION 


Consider  an  arbitrary  transmission  line  with  arbitrary  loads  at  both  ends  (see 
Fig.  1).  The  differential  equations  relating  the  voltage  V  and  the  current  I  along 
the  line  are  expressed  by 


.  jjj  —  c  +  G  V 
_  C0_+G0V 


where  L„  .  C„,  R„  and  G„  are  the  inductance,  capacitance,  resistance,  and  conductance 
per  unit  length,  respectively,  and  must  be  regarded  as  functions  of  frequency  to 
preserve  generality.  The  solutions  for  time-harmonic  excitation  are  usually  writ¬ 
ten  in  the  frequency  domain  as  (  «  =  2irf  is  angular  frequency ) 


—fix  fix 

v(aMt)  =  Ae  +  Be 


K«*x)  =  ±~  |ac  *x  -  Be**x  | 


where 


0  =  +  Y*U>  KG0  +  jwC0  ) 


+  j«l-0 


G„  +  jaiC0 


0  and  ZQ  are  complex  leading  to  an  attenuation  of  the  signal  as  it  propagates 
through  the  medium.  If  the  terminations  are  linear.  ( i.e. .  Zs(t)  and  ZL(t)  constant 
with  time),  the  coefficients  A  and  B  can  be  determined  by  matching  boundary  con¬ 
ditions  at  x=0  and  x=l;  next,  an  inverse  Fourier  transform  approach  can  be  used  to 
solve  for  the  tune-domain  solution.  On  the  other  hand,  il  the  terminations  are 
nonlinear  or  tune-changing,  then  the  boundary  conditions  must  be  formulated  in 
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the  time  domain  as 

vk(t)  =  V(t.O)+Z,(t)  I(t,0) 

(4) 

Zt<t)I<U)  =  V<U> 

(5) 

where  r,(t)  is  the  source  voltage.  Z*(t)  and  ZL(t)  indicate  the  time  variations  of  the 
source  and  load  terminations,  respectively.  For  any  time  greater  than  t.  Z,(t)  and 
ZL(t)  are  not  known  since  they  depend  on  the  voltage  and  current  solutions  at  time 
t.  In  tact,  the  evaluation  ol  Z,(t)  and  ZL(t)  may  involve  several  iterations  consist¬ 
ing  of  solving  the  propagation  equations  with  trial  values  until  convergence  to  the 
true  impedance  values.  Transforming  conditions  (4)  and  (5)  into  the  frequency 
domain  is  inappropriate  and  a  time-domain  formulation  is  thus  necessary.  Like¬ 
wise.  (2)  and  (3)  cannot  he  conveniently  inverted  into  the  time  domain  in  a 
closed-form  and  constrained  to  satisfy  (4)  and  (5).  This  limitation  arises  not  only 
because  of  the  dispersive  and  f requency -dependent  characteristics  of  the  line,  but 
also  because  the  evaluation  of  the  coefficients  A  and  B  requires  an  a  priori 
knowledge  of  the  time-dependent  load  functions.  A  formulation  in  which  the 
causality  of  the  boundary  conditions  is  implied  thus  becomes  necessary.  The  use 
ol  scattering  parameters  allows  one  to  dehne  the  properties  of  the  transmission  line 
independently  from  those  of  the  terminations;  consequently,  by  properly  combin¬ 
ing  load  and  line  relations,  a  closed-form  algorithm  for  the  solutions  can  be 


derived. 


III.  SCATTERING  PARAMETER  FORMULATION 


Any  linear  two-port  network  can  be  described  as  a  set  of  scattering  parame¬ 
ters  or  (S-parameters)  which  relate  incident  and  reflected  voltage  waves.  These 
waves  are  variables  which  depend  on  the  total  voltages  and  currents  at  the  two- 
port.  If  ideal  (lossless  and  dispersion  less)  transmission  lines  of  known  characteris¬ 
tic  (reference)  impedance  Znt  are  connected  to  both  ports  of  a  linear  network,  then 
the  voltage  waves  on  the  reference  lines  (see  Fig.  2)  at ,  b, ,  a2  ,  bj  are  defined  as  the 
incident  and  reflected  waves  from  port  1  and  port  2.  respectively.  The  scattering 
parameters  are  then  known  to  satisfy  the  frequency  domain  relation. 

faj  —  Sjj&i  S|2&2  (b) 

b2  “  SjiAi  +  S22&2  (7) 

S„  and  S22  are  regarded  as  reflection  coefficients  whereas  S,2  and  S2,  are  the 
transmission  scattering  coefficients  of  the  network.  The  total  voltages  in  ports  1 
and  2  are  respectively  given  by 

V,  =  a,+bj  (8) 

V2=a2+b2  (9) 

and  the  expressions  for  the  currents  are 


ai  _  b( 
3  2  b: 

Zrf(  Zr,f 


( 10) 

(11) 


A  relation  beiween  the  propagation  characteristics  of  a  transmission  line  and  the 
associated  scattering  parameters  can  then  be  easily  derived.  For  a  single  mode  of 
propagation,  it  can  be  shown  that  tor  a  given  transmission  line 
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(I  -p2)® 

l-p-v 


S|2  ~  *21 


(I  -ct)p 
l-p2**2 


(12) 


a  —  e 


„  Zo(«)-Z„f 

Z0(«)+Zref 


(13) 


where  «■(•)  and  Zo(«)  account  for  the  dispersive  and  lossy  behaviors  of  the  line. 
The  scattering  parameters  of  a  transmission  line  depend  only  on  its  electrical 
characteristics  and  are  not  influenced  by  the  source  and  load  voltages  at  the  termi¬ 
nations;  however,  the  overall  response  of  the  system  is  a  combination  of  line  and 
termination  responses  and  can  be  obtained  by  cascading  the  various  sections  of  the 
network.  We  can  then  write  the  time-domain  equations  relating  the  voltage  waves 
of  an  arbitrary  line  terminated  with  nonlinear  loads  (Pig  3).  We  get  (subscripts  1 
and  2  refer  to  near  end  and  far  end.  respectively) 


b,(t)  =  SuU^a.U)  +  S12(t)*a2<t)  (14) 

b,(t)  =  +  S22(t)*a2(t)  (15) 

where  •  indicates  a  convolution  in  the  time-domain.  The  scattering  parameters 
S,,(t) ,  S,,(t) ,  S2!(t) ,  S22(t)  are  the  inverse  transforms  of  the  frequency  domain  S- 
parameters  and  can  be  viewed  as  Green's  functions  associated  with  the  lime  domain 
response  of  the  transmission  line,  due  to  a  single  frequency  source  at  the  termina¬ 
tions.  I  he  load  conditions  at  the  near  and  far  ends  are  now  directly  expressed  in 
the  lime-domain  by  looking  at  the  flow-graph  representation  of  the  system  (see 
Fig-  3) 

a,(t)  =  r,(t)b,(t)  +  T,(t)g,(t)  (lb) 


,  a2(t)  =  T,(t)b,(t)  +  T,(t)g,(t) 


k/-/.  /s'"..'-/ ..'V'. ■.  •  •  . 


(  17) 


Figure  3.  Time-domain  circuit  and  flow-graph  representation 
of  arbitrary  transmission  line  (TL)  with  nonlinear 
terminations  at  time  t.  The  sign  *  indicates  a 
convolution  between  time-domain  scattering  parameters 
and  the  independent  voltage  waves  a}  and  a2  as  per 
Equations  (14)  and  (15)  . 


in  which  r,(t) ,  r2(t) ,  T^t) ,  T2(t)  are  the  reflection  and  transmission  coefficients 
associated  with  near  and  far  ends  respectively. 


Tj(t)  = 


Zi(t)  +  Zref 


i\(t)  = 


zs(t)  -  zrrf 
z,<t)  +  znr 


(18) 


In  Eqs.  (14)  and  (15),  each  of  the  convolutions  terms  can  be  expressed  as 

t 

V*}58  f  S|j(t— r)aj(r)dr  (19) 

o 

Since  the  algorithm  to  be  derived  must  be  amenable  to  computer  usage,  it  is  desir¬ 
able  to  discretize  ( 19)  as 


Sjj(t)  =  £  Sjj(t— r)aj(r)  At 

T~l 


(20) 


t-i 


S,,(t)  =  Si/0)  «j(t)Ar  +  £  Sjj(t — t)  aj(r)Ar 

r— 1 


or 


S,j(t)  =  S'i/OJa/t)  +  (21) 

i-i 

where  At  is  the  time  step  and  S'i,(0)  =  S,j(0)Ar.  H;j(t)  =  £  SSj(t — r)  aj(r)Ar 

r  =  l 

represents  the  history  of  the  line  and  depends  on  information  up  to  time  t-1. 
Causality  insures  that  the  aj  's  are  known  for  r<t,  which  allows  the  use  of  this 
information  tor  the  determination  of  the  aj's  at  r=t.  We  hrst  substitute  (21  )  into 
(  14 )  and  (  15 )  and  obtain 

b,(t)  =  S'u(0)  a,(t)  +  S',,(0)a:(t)  +  (22) 
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bj(t)  =  S'2I(0)  »,(t)  +  S‘22<0)  »2<t)  +  H2t(t)  +  H22(t)  (23) 

Combining  the  above  equations  with  those  for  the  forward  waves  (16)  and  ( 17), 
one  gets 


[l  -  r2(t«22(0)]  +  r,(t)M,(t)] 

l  -4 ) 

r^t^ufoKT^tig^t)  +  r>(t)M2(t)] 

+ - ^ - 

, .  [i  -  rjLtrstjm  [t^U)  +  r2(t)M2(t>]  . . . , 

• - m -  ( 25  ’ 

r^tJs’iitoXTiUJgift)  +  r,(t)M,(t)] 

+ - - 

A(t)  =  [i  -  r^tJSnfoilti  -  r2s22(o)j  -  r,(tis',2(o)r2(t)s'21«»  ( 26) 

where  M,(t)  =  Hu(t)  +  H,2(t)  ancjz  M2(t)  =  H>,(t)  +  H22(t).  b,(t)  and  b^it)  are 

recovered  using  (22)  and  (23).  and  the  total  voltages  at  ports  1  and  2.  by  using  (8) 
and  (9).  Figure  4  shows  the  flow  graph  representation  for  the  transmission  line  at 
time  t  in  which  the  memory  of  the  network  has  been  included  in  the  terms  M](t) 
and  M2(t).  The  independent  terms  are  g,<t)  and  g2(t).  M,(t)  and  M2(t)  are  also 
independent  and  contain  the  information  pertinent  to  the  history  of  the  line. 
Numerical  considerations  are  of  importance  for  the  simulation  since  they  determine 
the  speed  of  the  computations  involved.  Expressions  ( 24 )-( 26 )  can  be  further 
reduced  by  observing  that  for  transmission  lines  with  nonzero  length.  S'12(0)  and 
S',,<0)  must  vanish  since  a  hmte  duration  is  requirt'd  for  an  arbitrary  signal  to  pro¬ 
pagate  through  the  line.  The  above  relations  then  become 


ft 


:y,:. 


g,(t)  t,  (t) 


Sato 


M2(t) 


Mi(t) 


S,£o) 


02(t)  Tz(t)  g2(t) 


Figure  A.  Equivalent  flow-graph  representation  of  arbitrary 
transmission  line  with  nonlinear  time-varying 
terminations  at  time  t.  The  history  of  the  network 
is  Included  via  the  memory  variables  Mi(t)  and  M2(t). 
Note  that  this  representation  involves  only 
multiplications  In  the  time  domain. 


[1  -  r2(t»'22(0)][TI(t)g1(t)  ♦  y,(t)] 


[I  -  iVtfc  ^OMTjUWt)  +  y2(t)j 


bj(t)  =  SutOk^t)  +  M,(t) 


(27c) 


bz(t)  *  S  22(0)»2(t)  +  M2(t) 


(27d) 


Ait)  —  [1  -  rjtxs.^o)]!!  -  r2(tW22(o)] 


Computational  limitations  in  these  expressions  are  essentially  determined  by 
M,(t)  and  M2(t)  which  contain  the  history  of  the  network  and  involve  the  voltage 
wave  solutions  from  previous  time  steps,  in  the  case  where  losses  and  dispersion 
are  neglected,  the  frequency  domain  scattering  parameters  associated  with  the 
transmission  line  become 


S„(o.)  =  S22(«)=  (1— e--V-lp 

-2^ 

1  -Ve  v 


c<i  cii  ( 1  —  p- )  e  v 

S,2(m)  =  S21(«)  =  - - - ~ 

-2>~ 

1  ~  (St  v 


where 


Z„  Zr,f 

"h  Zrr, 


Since  Z„  the  c harat tens! it  impedance  ot  the  line  is  constant  with  time,  and  since 


Zn(  the  reference  impedance  is  arbitrary,  one  can  choose  Z,^*Z0  which  leads  to 
p  =  0  and 


Su(«)  =  Sjjftt)  =  0 

(32) 

S„(m)  B  Sjiftt)  =  e  v 

(33) 

Therefore,  the  time  domain  Green's  functions  associated  with  the  scattering  param¬ 
eters  are 


S„(t)  =  S22<t)  =  0 

(34) 

S.,(t)  =  S21(t)  =  S(t-1) 

(35) 

V 

M,(t)  =  a2(t—l/v) 

(36) 

M2(t)  =  a,(t— 1/v) 

(37) 

A(t)  =  1 

(38) 

We  then  obtain  simple  expressions  for  the  forward  and  backward  waves 


a|(t)  =  T,(t)gi(t)  +  r,(t)a2(t— 1/v) 

(39) 

a2(t)  =  T2(t)g2(t)  +  r2(t)a,(t— 1/v) 

(40) 

b,(  t)  =  a2(t— 1/v) 

(41  ) 

b2(t)  =  a,(t— 1/v) 

(42) 

The  advantage  ot  the  above  expressions  lies  in  their  computational  efficiency 
xinu'  only  a  search  is  involved  in  the  evaluation  of  the  history  ol  the  network  and 
no  summations  of  previously  calculated  terms  are  needed. 


IV.  MODELS  FOR  TERMINATIONS  AND  DEVICES 


Thus  far,  this  study  concentrated  on  simulating  the  time-domain  transient 
response  for  arbitrary  transmission  lines  terminated  with  nonlinear  devices.  In 
this  section,  we  examine  the  nature  of  the  terminations  and  the  manner  in  which 
they  are  to  be  represented  in  a  form  consistent  with  the  relations  derived.  Several 
techniques  are  available  that  convert  reactive  elements  and  nonlinear  devices  to 
time-varying  causal  resistances  as  well  as  voltage  or  current  sources.  We  briefly 
overview  two  of  these  techniques,  namely,  the  Backward  Euler  scheme  and  the 
Newton-Raphson  (NR)  algorithm.  A  more  detailed  development  can  be  found  in 
[14].  Since  formulation  and  solution  are  in  the  time  domain,  every  element  must 
have  an  equivalent  network  in  the  time  domain.  The  Backward  Euler  algorithm  is 
a  numerical  integration  scheme  and  its  use  in  representing  reactive  elements  in  the 
time  domain  is  illustrated  below.  Consider  a  capacitor  C.  with  a  current  voltage 
relation  given  by 


I 


(43) 


If  we  disc  ret  i/e  the  time  variable  by  choosing  a  time  step  h.  then  the  voltage  Vn+1 
at  tune  tn+1  =  (n+l)h  can  be  approximated  in  terms  of  variables  at  t  =  nh  as 

Vn+,  =  V„  +  hVn+l  (44) 

where  the  superscript  '  indicates  a  derivative  with  respect  to  time.  Making  use  of 
( 43  ).  we  get 


(45) 


Equation  (45)  tun  be  represented  bv  the  equivalent  linear  one-port  model  at  time 
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C  C 

tn+1  =(n+l)h  (see  Fig.  5)  with  conductance  G  =  —  and  current  source  J„  =  — Vn. 

n  h 

Such  transformations  can  be  applied  to  inductances  as  well.  The  efficiency  of  the 
method  depends  on  a  proper  choice  of  the  time  step  which  determines  the  stability 
of  the  numerical  solution.  Other  techniques  of  numerical  integration  can  be  used  if 
greater  stability  is  desired  [14]. 

Nonlinear  elements  such  as  diodes  and  transistors  must  be  reduced  to 
equivalent  networks  with  only  linear  elements  at  time  t.  If  the  nonlinear  current 
voltage  relation  for  these  elements  is  known,  then  an  iterative  scheme  such  as  the 
Newton-Raphson  algorithm  can  be  used  to  seek  a  solution.  The  circuit  representa¬ 
tion  of  the  Newton-Raphson  technique  is  illustrated  in  Figure  6.  At  a  particular 
time,  a  guess  value  for  the  voltage  is  chosen  to  which  a  current  is  associated  via  the 
I-V  relations  that  determine  the  operating  point  Q.  The  next  guess  is  then  relaied 
to  the  previous  one  by 

-i 

'in  “I*  —  -gj-  Ik  (46) 

At  each  iteration  step,  the  resulting  equivalent  circuit  is  composed  of  a  linear 
conductance  of  value  Gk  =  dl/dV  at  Vk  and  a  current  source  with  value  given  by 
Jk  =  I(Vk)— GkVk.  Solving  the  combined  transmission  line.  Newton-Raphson 
equivalent  circuit  problem  at  each  iteration  step  will  lead  to  the  actual  representa- 
Mon  ol  the  termination  at  time  t. 

Once  linearization  and  discretization  are  performed,  time  domain  values  are 
available  for  the  terminations  equivalent  impeda net's  or  generators.  These  expres¬ 
sions  are  causal,  since  their  values  at  anv  time  t  depend  on  the  history  of  the  net¬ 
work  which  renders  impossible  an  ■  i  priori  knowledge  of  the  tune  variations  ol  the 
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V.  APPLICATIONS  AND  PRACTICAL  CONSIDERATIONS 

Many  applications  in  microwave  and  digital  communications  require  the  use 
of  transmission  lines  terminated  with  nonlinear  devices.  Distortion  and  noise  arise 
when  the  terminations  are  not  matched  to  the  line  impedance.  Moreover,  if  losses 
and  dispersion  are  present  in  the  line,  attenuation  and  delay  come  into  account. 
The  combination  of  these  effects  needs  to  be  modeled  and  simulated  on  a  reliable 
CAD  tool.  The  above  algorithms  were  used  to  predict  signals  at  the  near  and  far 
ends  of  a  73  ft  stripline  structure  terminated  with  ALS  inverters.  The  length  of 
the  line  was  50  inches.  The  pulse  characteristics  are  summarized  in  Figure  7. 
Results  were  compared  with  experimental  measurements  and  good  agreement  indi¬ 
cated  the  validity  of  the  assumptions. 

Equations  (27)  were  used  to  predict  the  waveforms  on  a  microstrip  line  hav¬ 
ing  dispersive  and  lossy  characteristics.  Predicted  distortion  waveforms,  line 
parameters  and  simulation  data  are  summarized  in  Figure  8.  The  dispersion  model 
of  Bhartia  and  Pramanick  [3]  was  used.  This  model  provided  the  frequency 
domain  data  with  a  frequency -dependent  effective  dielectric  constant  and  charac¬ 
teristic  impedance.  A  Vf  behavior  was  assumed  for  the  resistance  per  unit  length  to 
model  the  skin  effect.  Waveform  attenuation  as  well  as  rise  and  fall  time  degrada¬ 
tion  can  be  observed  from  the  simulations.  Wafer  scale  integration  and  VLSI  inter¬ 
connections  that  possess  higher  loss  parameters  [15]  should  be  expected  to  present  a 
higher  level  of  signal  degradation. 
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Figure  7.  Comparison  of  theoretical  (plots)  and  experimental 
(photographs)  simulations  for  a  strlpllne  structure 
terminated  with  ALS  Inverters  at  near  (left)  and 
far  (right)  ends.  The  line  length  is  50  Inches. 

The  low-frequency  characteristics  are  Z0  -  73ft  and 
v  ■  0.169/ns.  Pulse  width  is  160  ns;  rise  and  fall 
times  are  4  ns. 
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VI.  CONCLUSION 


This  study  explored  some  important  aspects  of  interconnections  for  digital 
and  microwave  applications.  The  problems  of  losses,  dispersion  and  load  non- 
linearities  were  analyzed.  A  closed-form  algorithm  was  derived  for  the  simulation 
of  an  arbitrary  time-domain  signal  on  a  structure  having  all  the  above  properties. 
The  algorithm  assumed  that  the  frequency-dependent  characteristics  of  the  line 
were  available  as  well  as  large-signal  models  for  the  terminations.  Future  work 
includes  the  derivation  of  a  suitable  loss  and  dispersion  model  and  the  extension  of 
the  algorithm  for  modeling  n-line  multiconductor  systems. 
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